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ABSTRACT 

We present high spatial and spectral resolution observations of sixteen Galac- 
tic compact and ultracompact H II regions in the [Ne II] 12.8 /im fine structure 
line. The small thermal width of the neon line and the high dynamic range of 
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the maps provide an unprecedented view of the kinematics of compact and ul- 
tracompact H II regions. These observations solidify an emerging picture of the 
structure of ultracompact H II regions suggested in our earlier studies of G29.96 
and Mon R2; systematic surface flows, rather than turbulence or bulk expansion, 
dominate the gas motions in the H II regions. The observations show that almost 
all of the sources have significant (5-20 km s~^) velocity gradients and that most 
of the sources are limb-brightened. In many cases, the velocity pattern implies 
tangential flow along a dense shell of ionized gas. None of the observed sources 
clearly fits into the categories of filled expanding spheres, expanding shells, filled 
blister flows, or cometary H II regions formed by rapidly moving stars. Instead, 
the kinematics and morphologies of most of the sources lead to a picture of H II 
regions confined to the edges of cavities created by stellar wind ram pressure 
and flowing along the cavity surfaces. In sources where the radio continuum and 
[Ne II] morphologies agree, the majority of the ionic emission is blue-shifted rel- 
ative to nearby molecular gas. This is consistent with sources lying on the near 
side of their natal clouds being less affected by extinction and with gas motions 
being predominantly outward, as is expected for pressure-driven flows. 

Subject headings: interstellar medium - ultracompact H II regions - fine structure 
line 



1. Introduction 



The short Kelvin-Helmholtz timescales of massive young stellar objects (YSOs) allow 
these objects to begin emitting ionizing ph otons and to excite nearby material while still 
embedded in their natal clouds (lLadalll999l ). One of the earliest observable manifestations 
of newly formed OB stars is the ultracompact HII (UCHII) region. These embedded HII 
regions are initially small (D< 0.1 pc) and dense (?7,e~10^ cm'^) and have a typical emission 



measure (EM) ~ 10^ pc cm ^ fjChurchwelll |2002| ) . As HII regions evolve their sizes increase 



but they are still considered to be UCHII regions as long as their emission measures (or 
surface brightnesses) remain high. 

Since most stars with M > 8 M© should form UCHII regions shortly after their birth, 
studies of UCHII regions can help us learn about physical conditions in high-mass star form- 
ing regions and the evolution of the physical conditions as the massive stars turn on. A large 
percentage, > 50%, of UCHII regions have cometary, core-halo or irregular morph ologies 
dWood fc ChurchwelllE989bl : lKurtz e^ll994J : Iwalsh et allEggsilde Free et al.ll2005l) argu- 
ing against the classical expanding sphere model of HII regions (ISavedoff &: Greenelll955l ). 
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In addition, statistical arguments based on the Galactic population of UCHII regions lead 
to estimates of the lifeti nae of these regions that are much longer tha n expected from the 
expanding sphere model (jWood &: Churchwelllll989al : iHoare et al.ll2007l'). Man y models have 
been proposed for UCHII region formation (see the review by lChurchwelllll999l . and references 
therein), with solving the lifetime problem a major goal. However, the age of an individual 
UCHII region is hard to determine through observations. A number of theoretical models 
have been geared toward explai ning both these long apparent lifetimes and the observed 



morphologies of UCHII regions (IHoUenbach et al.l 1 19941 : iDyson et al.l Il995l : iRedman et al. 
19961 : Iwilliams et al.lll996l : iRedman et al.lll998h . 



Observations of the morphology and gas kinematics of UCHII regions can help us to 
constrain models for the structure and evolution of the regions much more definitively than 
observations of the morphology alone. The morphological type is not a unique i ndicator of the 
UCHII region physics because many models can produce similar morphologies (ITenorio-Tagle 
1979 1: lYorke et al.1ll983l : iMac Low et~allll99ll : Ivan Buren fc Mac Lowlll992l : krthur fc Hoare 
20061 ). While it is difficult to choose between different models based on the morphologies 
of HII regions, their kinematic signatures are very distinct. To understand the structure 
and evolution of the regions around the youngest massive stars, it is therefore necessary to 
observe gas kinematics inside UCHII regions with the highest spectral and spatial resolu- 
tion. The large extinction associated with molecular clouds around UCHII regions, however, 
prevents direct obs ervation at optical wavelengths and even at near-infrared wavelengths 
(JDopita et al.ll2006l ). Thus, observations at longer wavelengths become crucial to UCHII 
region studies. 

Radio interferometric observations can provide data with high angular and spectral reso- 
lution. At high enough frequencies, where non-LTE effects become less important, hydrogen 
radio recombination line (RRL) emission has the same dependence on the electron density 
as thermal free-free radio continuum emission, and the emissivities of these lines are well 
known. RRL observations can provide information on gas density and kinematics, and, to 
a limited degree, on gas temperature. However, thermal broadening of hydrogen RRLs due 
to the low mass of hydrogen atoms degrades the spectral resolution of radio spectroscopy. 
The thermal line width of the RRLs (or any hydrogen recombination lines like Bra and 
Br7, for that matter) in an ionized region (Te~10^ K) is ~ 20 km s~^, comparable to the 
velocities of bulk flows in many sources. Furthermore, the incomplete uv sampling of radio 
interferometric survey can both overemphasize compact structures and, more generally limit 
the spatial dynamic range of the observations. 

Mid-IR fine structure line observations offer an alternative to radio and IR recombina- 
tion line observations. Fine structure lines have small thermal line widths because they are 
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emitted by heavy elements. The importance of this feature of the fine structure hues has 
been illustrated previously by high spectral resolution maps of [Ne II] emission in Mon R2 
and G29.96-0.02 which show single lines or multiple components of lines with widths sig- 



nificantly smaller tha. n the thermal widths of the hydrogen recombination lines (jJaffe et al. 



20031 : IZhu et al.ll2005l ). In some cases, the broad, largely Gaussian hydrogen recombination 
lines are, in fac t, superposi t ions o f several much narrower kinematic components (see espe- 
cially Fig. 3 of IJaffe et al.l (120031 )). In addition, although the emitting species are several 
orders of magnitude less abundant than hydrogen, the fine structure lines are brighter than 
recombination lines because collisional excitation cross-sections are much larger than the 
relevant recombination cross-sections. Typical electron densities over all compact and most 
ultracompact H II regions are lower than the critical density of the [Ne II] 12.8 /im transi- 



tion (^5 X 10^ 



cm 



so this line shares the property of the recombination lines of having 



a line strength proportional to emission measure. We should mention that turbulence can 
also have significant effects on line widths. Moreover, turbulence affects both RRL and fine 
structure line observations and prevents detections of velocity variations (due to gas bulk 
motions) with scales smaller than the intrinsic line width combining thermal and turbulent 
effects along any line of sight direction. 

Our previous [Ne II] line observations of a few UCHII regions have demonstrated the 
effectiveness of usin g high spectral and spatial resol ution maps of fine structure lines to 
study such regions ( jJaffe et al.ll2003l : IZhu et al.ll2005l ). In these studies, we demonstrated 
that the kinematic patterns in the ionized gas in Mon R2 and G29.96-0.02 are consistent 
with tangential flows along the surfaces of thin, more-or-less parabolic shells. The goal 
of the current paper is to use [Ne II] mapping as a tool to examine a larger sample of 
compact and ultracompact HII regions with a range of shapes and sizes in order to look for 
a common thread in the kinematics that might lead to a better understanding of the physics 
and evolution of these sources. We present scan-mapped spectroscopic observations of fifteen 
Galactic compact and ultracompact HII regions. This survey is the first mid-infrared high- 
spectral resolution study of gas kinematics in a large sample of UCHII regions. 

In § m we describe our observational method. We introduce our data reduction pro- 
cedure in § [51 In § 31 we investigate the morphology and kinematics of individual objects. 
Finally, the conclusions drawn from our observations are presented in § O 



Observations 



Candidates in our surve y were selected based on their radio continuum flux density levels 
( IWood &: Churchwelllll989bl ). biasing our list toward radio-bright compact and ultracompact 
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HII regions. Observations were performed with the high spectral resolution, mid-IR (5- 
25 ^ni) cross dis persed spectrograph TEXES (Texas Echelon Cross Echelle Spectrograph, 



Lacy et al.ll2002l ) on the 3.0 meter NASA IRTF on Mauna Kea in Hawaii. 



We list our observational parameters in Table [T] Our spatial and spectral resolution 
were approximately 1.8" and 4 km s~^. For each object, a north-south oriented slit was 
used to scan the object from west to east. We chose the step size along the scans so that 
there were at least two steps per slit width and chose the scan lengths to cover the emission 
seen in radio continuum maps available in the literature. Before each scan, we offset the 
telescope to the west several arcseconds beyond the edge of the maps shown in the figures. 
Several extra steps were taken at the beginning and the end of each scan in order to measure 
the sky background. The lengths of the scans and the initi al offsets were ad justed after 
examining a preliminary scan with our quick-look program ( Lacy et al.l |2002| ) . We made 



multiple overlapping scans, offset north-south, to cover the sources. TEXES takes ambient 
temperature blackbody and sky flats at the beginning of each scan for flux calibration. 



Data Reduction 



We use a custom Fortran reduction program (iLacy et al.l |2002| ) to correct optical dis- 
tortions, fiat field, remove cosmic ray spikes, and fix bad pixels. The same program also 
does wavelength and flux calibrations. We use atmospheric lines for wavelength calibration. 



with wavelengths obtained from HITRAN (JRothman et al.lll998l ). We use an IDL script to 



subtract sky background and combine multiple scans. The sky emission at each slit position 
is linearly interpolated from the sky frames at the beginning and the end of each scan. Some- 
times, we start scans too late or end scans too early. Extended line emission is present in 
our sky frames and is subtracted off from other frames. This will result in an underestimate 
of line flux from targets. Multiple scans are cross-correlated, shifted and added to create 
the datacubes. Normally, no shift is needed along the spectral direction. One exception was 
the source W3, whose datacube was obtained by combining data from two nights and where 



the m aps were registered using the compact continuum source IRS5 (jWynn- Williams et al. 



19721 ). Both the Earth's motion and the spectral settings of the instrument on two different 
nights affect the wavelength calibration. The continuum flux level at each pixel in the map 
is obtained from spectral points off the line and is subtracted from the spectra. Finally, we 
resample the datacubes so that each pixel has a size of 0.3"x0.3" on the sky. The coordinates 
of the resulting continuum maps and [Ne II] line maps are relative to the peak positions in 
the integrated [Ne II] line maps. When radio continuum maps are available, they are cross- 
correlated with the [Ne II] maps so that both maps can be plotted on the same scale. The 
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coordinates in Table 1 were taken from the literature and served as reference points during 
our line observations. 



4. Results 

In this section we discuss highlights of individual objects and features common to several 
sources. [Ne II] line emission datacubes for all objects are available online. Note: the data 
cubes will be available through the ApJ web site. 



4.1. G29.96 -0.02 

G29.96 -0.02 (G29.96 hereafter) is a cometary UCHII region at a distance of ~ 7.4 kpc, 
is ~20"xl5" in size, and is one of the most luminous UCHII regions in the sky at 2 cm 
wavelength. Th e ionized region is co nsidered as an archetype of the UCHII regions caused 
by bow shocks (JMac Low et al.lll99ll ). although some past observations of the gas velocities 
along the symmetry ax is of this source argue ag ainst this interpretation, implying that the 
source is a blister flow flLumsden fc Hoard 1 19961 ) or a blister flow modified by a stellar wind 
(JLumsden fc Hoard 1 19991 ). We presented [Ne II] line m apping observati ons and compared 
gas kinematics in the region with a bow shock model in IZhu et al.l (120051 ). We include here 
a summary of the model and results for G29.96 from the earlier paper as a template for 
understanding the kinematics of some of the other sources in the sample. The bow shock 
model predicts that an approximately paraboloidal shell of swept-up material will form in 
front of a star moving supersonically inside a molecular cloud. The shell is supported by 
the ram pressure of the stellar wind and the ambient medium. Swept-up material, including 
ionized and neutral gas, is forced to flow along the surface of the shell. With different 
viewing angles for the shell, the model can reproduce many of the observed features in ionic 
line observations of G29.96 as well as several other UCHII regions in our sample, including 
morphologies and the distinctive shapes that parabolic flows have in position-velocity (p-v) 
diagrams. 



Figure [T] shows a 2 cm radio continuum map from iFey et al.l (119951 ) and our spectrally 
integrated [Ne II] line map of G29.96. The two maps closely resemble each other, indicating 
that most of the [Ne II] emission arises from gas with n<nr.r if. and that the e xtinct ion to the 
gas is either rather uniform or small at 12.8 /im. In fact, JMorisset et al.l ((2002) observed 
that the electron density in th e region is sub-critical (Ne~680 cm~^ - 5.7x10^ cm~^) and 
Martin-Hernandez et al.l (12002! ) reported low infrared extinction (A^^^l.G) toward the region. 
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Our findings are consistent with these observations. Both maps show an emission arc open 
to the east with extended emission in the eastern part of the region. A filament extends from 
the southern edge of the nebula to the east, forming a "tail" of the cometary HII region. 
Channel maps at nine equally spaced velocities are shown in Figure [21 These channel maps 
show a morphology change from a thin crescent in the most blue-shifted channels to a fan- 
like region in red-shifted channels. Figures [3] and H] show position-velocity (p-v) diagrams 
of several representative cross cuts made through the observed and model data cubes. The 
observations and the model display similar p-v structures (as shown in our previous paper, a 
characteristic "A" shape for cuts perpendicular to the nebular symmetry axis and a "7" shape 
for cuts parallel to the axis). Similar changes in structure from cut to cut can also be seen. 
The similarities in line emission morphology and gas kinematics between the observations 
and the model indicate that the ionized gas in G29.96 flows along a thin shell at ~ 20 km s~^, 
as it would if it were produced by a bow shock. We use these distinctive patterns in the p-v 
diagrams as indicators of the presence of surface flows in other sources where we have not 
carried out detailed model fits. However, one clear discrepancy between the bow shock model 
that gives the correct shape to the contours in the model p-v diagrams in Figures [3] and H] 
and the observed sources is the ~15 km s~^ difference between the observed and modeled 
velocity of the ambient neutral material. We discuss in § [5] modifications of the bow-shock 
model to explain this discrepancy. 



4.2. G5.89 -0.39 



G5.89 -0.39 a lso called W28 A2(l), is 2 .6 kpc from the Sun flDownes et all Il980 



Kim fc Kodl2003l ). IWood fc Churchwelll (Il989bl ) measured continuum images of the source 
at both 2 cm (Figure E], left) and 6 cm. A shell-like morphology with significantly limb- 
brightened walls surrounding the slightly elongated central cavity can be recognized. The 
G5.8 9 region also hosts a source of extremely broad (FWHM ~140 km s~^) CO J =l-0 emis- 
sion (jShepherd fc Churchwelllll996l : lHarvey fc Forveillelll988l ). lAcord et al.l (119981 ) have used 
multi-epoch radio continuum observations to argue tha t the ionized shell is expand ing at 
~35 km s~^. Radio recombination line observations of iRodriguez-Rico et al.l (120021 ) show 
a complex velocity structure across the region, which cannot be explained with a consis- 
tent dynamical model. Mid-IR dust continuum observations show that G5.89 deviates from 
spherical symmetry and resembles a shell partially open to the south, with significant SW 
and SE protrusions from the geometric center, implying an anisotropic expans ion of the 



ioniz ation front due to density variations in the surrounding molecular medium (JBall et al. 
1992h . 



Our [Ne II] 12.8 yum line map (Figure [5], right) shows a morphology similar to that of 
the mid-infrared continuum, with a bright arc along the north side of the source. With our 
~ 1.5" resolution, the kinematic signature of an expanding shell is not apparent, although 
there is some evidence that the line is broadest toward the center of the source (Figure [H]). 
The majority of the emiss ion is blue-shifted by ~9 km s~^ with respect to the molecular 
cloud (JHatchell et al.l Il998l ) , and the extent of the region becomes relatively smaller toward 
the red end of the spectrum. 



4.3. Gil. 94 -0.62 



Located at a distance of 4. 2 kpc (IChurchwell et al.lllQQOl ). Gil. 94 -0.62 was categorized 
by IWood fc Churchwelll (jl989bl ) as a cometary UCHII region. In their radio continuum map 
(Figure[7l left), the region shows an emission arc with its apex pointing to the west and diffuse 
emission ex tending to the east and north . A broad ^^CO(J=1-0) line was observed toward 
the object (IShepherd fc Churchwelll 119961 ). The CO line width (full width at zero int e nsity) 
of over 28 km s~^ suggests the presence of high velocity outflows. iDe Buizer et al.l (120031 ) 
imaged the region at 11.7 /im and 20.8 yum. The region does not show a cometary shape in 
their mid-IR images and, although their pointing has an astrometric accuracy about 1", the 
radio peak does not match any structure in the mid-IR images. The overall appearance of 
the region at mid-IR and radio continuum wavelengths does not match well. 

O ur [Ne II] 12.8 um line map (Figured right) matches the mid-IR continuum observa- 
tions (|De Buizer et al.ll2003l ) very well. By registering our map with the continuum map in 
De Buizer et al.l ( l2003l ). we determine the position of the radio continuum peak in our map, 
marked with an 'X' at (-1.17", 1.67"). Five separate emission peaks (designated as A-E) 
can be seen in ou r integrated [Ne II] map . They all have corresponding 11.7^m and 20.8^m 
emission peaks in lDe Buizer et al.l (120031 ). Linear sizes of these clumps range from ~2" in C 
to over 7" in A and E. Like the mid-IR continuum, the [Ne II] morphology is a poor match 
to that of the radio continuum. 

The spectra of the clumps C, D and E each can be fit with a single Gaus sian profile with a 



cente r velocity very close to the molecular cloud velocity (Vlsr —39 km s~^. IChurchwell et al. 
1990l ). Figure [8] shows line profiles in the region of A, B and C. Over this part of the source, 
the lines are broad, with most of the additional emission to the red of the cloud velocity. 
The profiles in Figure [8] show that most of the [Ne II] emission occurs over the range Vlsr 
= 30-60 km s~^. Complex morphology, probably induced by dust extinction, prevents a 
detailed analysis of the kinematics in the region. 
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4.4. G30.54 +0.02 

G30.54 +0.02 has a limb-bri ghtened, or "horseshoe" shape in both the radio continuum 
map (jWood fc Churchwel]||l989bl ) and our [Ne II] line map (Figure[9]). The more concentrated 
emission with a bigger opening in the radio continuum map probably is the result of the 
higher resolution of the interferometer observations and of the uneven sensitivity of the 
snapshot mode aperture synthesis to compact and extended emission. Both radio and [Ne II] 
line observations show that the emission arc breaks into at least three components. 

Figure [10] shows [Ne II] profiles across the region while Figure [11] displays position- 
velocity diagrams on selected cuts along and across the symmetry axis. Large line widths 
(20-30 km s~^ FWHM), together with smooth profiles imply that turbulence or substantial 
small-scale motions are present and may mask the bulk motion of the gas. Nevertheless, a 
change of velocity pattern within the structure is still apparent. 

The velocity pattern along the bright arc is consistent with significant tangential flows 
of the type seen in G29.96, but here mostly in the plane of the sky. At the vertex of the 
horseshoe, there is a single broad line with a velocity centroid close to the molecular cloud 
velocity (Figure [TT]) . On the arm of the horseshoe to the east of the vertex the line peak 
shifts ~ 7 km s~^ to the blue of the molecular cloud velocity while on the western arm the 
molecular line and [Ne II] peak velocities agree. [Ne II] emission is present over at least 
40-50 km s~^ at all positions (Figures [TO] and [TTj) . 

Along the symmetry axis and away from the bright arc, the match to a tangential flow 
is less clear. Rather than breaking up into two narrow lines, the emission on the axis and 
behind the vertex is at a single velocity close to that of the molecular cloud. A small velocity 
gradient across the apex may indicate an inclination to the plane of sky. 



4.5. G33.92 +0.11 

G33.92 +0.11 is particularly remarkable because of the apparent disconnect between its 
nondescript morphology and its dis tinctive position-velo c ity stru cture. It was first catego- 
rized as a core-halo UCHII re gion by] Wood fc Churchwell] (ll989b] ) and later as a shell UCHII 
region by ]Fey et al.] (1992]). ]Fey et al.l ( 19921 ) showed that the 20 cm extended emission is 
roughly shell-like, although the cavity in their 20 cm image is offset to the southeast of the 
emission peak, and the rest of the region has a morphology close to a core-halo region with 
cometary extended emission. C^^O observations s howed two emission co res, and one of them 
is associated with the peak of the UCHII region (jWatt fc Mundy|]l999l ). 
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The integrated [Ne II] line map shows an irregular peak surrounded by extended emission 
(Figure [T2I) . Mid-infrared continuum observations show similar structure to that seen in the 



[Ne II] line observations (JGiveon et al.l 120071 ). The major structure is a resolved ellipsoid 
oriented northeast-southwest with an additional extension to the east along the minor axis. 
The [Ne II] and radio continuum maps look very similar, but the central peak splits into two 
components in the neon line, possibly as a result of foreground extinction. 

Position-velocity diagrams along a few cross cuts through G33.92 (FiguredS]) show kine- 
matics identifiable with the kinematics typical of a flow along a surface similar to that seen in 
G29.96 including the characteristic "7" and "A" patterns, despite its more nondescript mor- 
phology. Profiles at individual positions are broad and frequently asymmetric (Figure [T4|) . 
Channel maps (Figure \T5\i show that three prominent components form an arc at blue ve- 
locities and persist to the ambient material velocity, Vlsr —108 km s~^. This arc has a 
sharp falloff to the west. A fourth emission peak and extended emission can be seen to the 
east of the arc. The positions of these emission peaks shift from channel to channel. This 
position shift is especially prominent for the peak to the east of the arc. It moves almost 
5" southeast from a position close to the arc in the channel at Vlsr—^^^ km s~^ to (+6,-3) 
in the channel at VisR—89 km s~^. The three components along the arc move a smaller 
distance to the northwest. Towards the red end of the spectrum, the emission components 
merge together and fade away. We note that the overall shape of the region changes from a 
thin crescent at blue-shifted velocities to a fan at red-shifted velocities. Bow shock models, 
like that for G29.96 predict a similar morphological change. 



4.6. G43.89 -0.78 



G43.89 -0.78 was classified as a cometary UCHII region by lWood fc Churchwelll (Il989bl ). 
Mac Low et al.l (Il99ll ) used G43.89 as one of their bow shock model tests. They concluded, 
based on the continuum maps that they nee ded a viewing angle between 4 5° and 60° from 
head-on in order to fit the morphology well. IShepherd &: Churchwelll ( 119961 ) found evidence 
for high- velocity molecular material in the region with a CO line width over 39 km s~^. Low- 
resolution 21 /im imaging observations showed the region as a point source with extended 
emission to the southeast (jCrowther fc Contill2003l ). while the radio continuum observations 
showed extended emission northwest of an emission peak. 



The cometary morphology and kinematics of the [Ne II] line emission (Figure [T6lfT7|) are 
consistent with flows along the surface of the region l ike those seen in G29.96. We fitted the 
p-v diagrams with a model like that used for G29.96 (IZhu et al.ll2005l ). The best fit was with 
a stellar motion of 15 km s~^ and 140° away from the direction to the observer. However, 
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as with G29.96, there is a discrepancy between the model, which predicts that the ionic 
emission should be redshifted relative to the molecular emission, and the observed emission, 
which is predominantly blueshifted. 



4.7. G45.07 +0.13 



G45.07 +0.13 was classified a s an UCHII region with a shell morphology by lTurner fc Matthews 

1984J). At a distance of ~6.0 kpc fichurchwell et alJEoOol ). it has a diameter ~0.9" at A=2 



cm 



( IWood fc Churchwelllll989bl ). The slight vel ocity gradient in t he H76a line was interpreted 



as the result of eith er an expanding ring (JGaray et al.l Il986l ) or ionized bipolar outflows 
( Lim fc White 19991). A group o f water masers lie l!' t o the north of the radio continuum 
peak jHofner fc Churchwel3ll99d : IPe Buizer e-TaDboOsh . 



G45.07 is unresolved in our 12.8 yum [Ne II] map at 2" spatial resolution (Figure [T5| 
right). The [Ne II] line toward G45.07 is m uch narrower (Vpy /if/vf ~ 20 km s~^) than is the 



H76a radio recombination line (~48 km s ^. iGaray et al.lll986l ). The difference in widths can 



arise either from pressure broadening of the H76a line or from the presence of high velocity 
material at very high density where the recombination line emissivity is proportional to n^ 
while the [Ne II] emissivity only scales with Ug. Dust continuum emission is observed at 
G45.07 and at a second source ~2" north of G45.07 (Figure [T51 left). The location of the 
second source matches the positions of the H2O masers. We did not detect any line emission 
from t he second source, supporting the suggestion that th is source is at an early evolutionary 
stage jHofner fc Churchwelll ^9961 : be Buizer et allboosh . 



4.8. G45.12 +0.13 



G45.12 +0.13 (G45. 12 hereafter) is located at a di stance of 6.9 kpc (jChurchwell et al. 



19901 ) . It was observed by IWood fc Churchwelll (jl989bl ) with the VL A at 2 and 6 cm and 
revealed at least 3 peaks with a roughly arc-like overall morphology. 

Our [Ne II] line observations do not completely resolve the three components along the 
radio emission arc, which is only ~3" long (referred to as Source "N" in Figure [T^ . The 
channel maps (Figure [201) show that the "horns" of the arc-like structure are more prominent 
at lower redshifts while the bright spot at the middle of the arc dominates at higher velocities. 
The p-v diagrams of a cut parallel to the symmetry axis of the source N do show "7" patterns 
of cometary H II regions (Figure [?H top panels), but a cut perpendicular to the axis does 
not show the expected "A" shaped pattern (Figure [211 bottom panels). It is not clear how 
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much the kinematics in G45.12N resembles that in a cometary HII region, because its "7" 
type p-v diagram is not dramatic compared to that in G29.96 (Figure [3]). 

We also find two fainter sources to the south of the arc which were not seen in the 
sparsely sampled VLA snapshot image. To show these two faint sources clearly, we extend 
contour levels down to 1% of the peak value in our line map. The southeastern source (SE) 
is extended and almost round-shaped. It is hard to determine its symmetry axes. Its p-v 
diagrams along RA and DEC cuts (Figure [22]) show evidence of gas bulk motion and most of 
emission is blue-shifted relative to the ambient molecular gas. Tentative "7" shaped emission 
distribution can be seen in these diagrams, suggesting that surface fiows may be present in 
the source. The southwestern source (SW) is irregular and contains multiple peaks. These 
emission peaks surround an emission minimum at the same location ~(2."7, -8") in almost all 
channels (Figure 1201) . A broken shell-like morphology with an opening to the west suggests 
an expanding shell structure for SW, but the ring of the line emission does not collapse at the 
ends of the Doppler shift-range as would be expected for an expanding shell. The velocity 
distributions in both N and SW are more or less symmetric about the ambient material 
velocity, which is 59 km s~^. Line emission from SE is relatively blue-shifted, and the width 
of the [Ne II] line is narrower in SE than in the other two sources. 



4.9. G45.45 +0.06 



G45.45 +0.06 was categorized by IWood fc Churchwelll (Il989bl ) as a cometary UCHII 
region. In the 6 cm radio continuum image, the source extends ~ 9" EW and 6" NS and 
has a sharp and bright ionization front to the north and extended emission to the south 
(Figure [231 left panel). Multiple emission peaks arranged along the ionization front form a 
distorted ring around a cavity. At a distance 6.6 kpc, t he integrated continuurn flux density 



i mplie s an 07.5 main sequence star as the ionizing star (jChurchwell et al.lll990l ). iFeldt et al. 



( 119981 ) acquired images at H and K' bands and three mid-infrared bands centered at 3.8^m, 
10.5y[xm and 11.7/um and found that the region contained a young OB cluster. 

The overall morphology of the [Ne II] line emission (Figure [231 right) does not look 
like that of the radio continuum map. Most of the line emission comes from an elongated 
structure aligned northeast to southwest w i th the extended emission to the southwest roughly 
matching that seen in Br7 by lFeldt et al.l (119981 ). The intermediate velocity (43-53 km s^^) 
channel maps (Figure [Ml) most nearly match the radio emission. The NE-SW elongated 
structure is primarily present in red-shifted velocity c hannels. Extende d emission north of 
the radio arc and a overlying compact 11.7 /xm source (JFeldt et al.lll998l ) dominate the blue- 
shifted channels. The differences between the [Ne II] and radio distributions could be due 
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to incomplete uv coverage of the snap shot radio observations or/and to greater sensitivity 
of [Ne II] line observations. 



4.10. W51 IRS2 



W51 IRS2 (W51d) is aluminous (L = 2-4 xlO^ L^. lErickson fc Tokunagalll980l : ljaffe et al 



19871 ) ultracompact HII region/molecular cloud core within a more extensive star forming 
complex (Carpenter & Sanders 19 98). Proper motion studies of HoO maser features plac e 
W51 IRS2 at a distance of ~ 7 kpc jOenzel et al1ll98ll : ISchneps et alj[l98ll : llmai et alJl2nn2l ). 
Elsewhere, we have pres ented 0.5^^ resolut ion maps of the [S IV] 10.5 /xm and [Ne II] 12.8 /im 
lines toward this source (JLacy et al.ll2007l ). These maps reveal a number of embedded UCHII 
regions, which dominate the [Ne II] emission, a more extensive surface flow, seen mostly in 
[S IV], and a prominent ~ 100 km s~^ jet, which we attribute to ionization of one lobe of 
a protostellar jet. The surface flow has a position-velocity pattern indicating that the gas 
is moving along a shell with its vertex pointed away from the observer. At the shell edges, 
the emission is blue shifted by 10-15 km s~^ with respect to the ambient cloud. The flow 
is less evident in [Ne II] where several compact sources are much brighter, either because of 
the lower extinction at 12.8 /im or because of the lower ionization potential of neon than of 
S++. 

We include our somewhat more extensive 2" resolution IRTF map of W51 IRS2 in the 
[Ne II] line here for completeness. Figure [25] compares the distribution of radio continuum 
and [Ne II] emission. The channel maps in Figure [26l as well as the top p-v diagram 
in Figure [57] show that the two components of the eastern source dl have velocities that 
differ by ~ 10 km s~^, with the southern component blueshifted by ~ 12 km s~^ from the 
molecular cloud velocity. The molecular emission spectra of CS show a second feature close 
to the [Ne II] velocity in the southern component (49.5 km s~^. Plume et al. 1997). The 
northern component of W51dl and W51d have [Ne II] velocities close to the velocity of the 
main component of the molecular cloud. We also notice that the ionized gas toward the 
center of the region is somewhat redshifted with respect to ambient material, which is likely 
because the ionized gas is pushed into the cloud by stellar/cluster winds. 



4.11. G61.48 +0.09B 



G61.48 -I-0.09B bel o ngs to the emi s sion n ebula complex Sh2-88B at a distance of 5.4 kpc 
( IChurchwell et al.lll990l ). I Gar ay et al.l (119941 ) compared kinematics from hydrogen recombi- 
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nation line and molecular line observations (JMacDonald et al.lll98ll : IChurchwell et al.lllQQOl ) 
and concluded that the region was undergoing a blister flow to the southwest. They argued 
that the hydrogen recombination line kinematics and the radio continuum morphology of 
the region were best explained by expansion of ionized gas in a non-uniform medium. Ra- 
dio continuum observations show that the region consists of two components, designated as 
Bl and B2 in Figure [281 The more compact and spherical B2 is located ~11" east of the 
peak of Bl, which is extended and has slightly curved contours that hint toward a cometary 
morphology. Our [Ne II] line emission distribu tion (Figure [22]), which is very similar to the 
distribution of Br7 emission (jPuga et al.ll2004l ). shows significantly different structure than 
the radio continuum map. The morphological difference between the [Ne II] line emission 
and the radio continuum emission is likely due to dust extinction. The [Ne II] emission has 
complex morphology and kinematics. Multiple components are present over a 35" by 20" 
area. A bright ridge of emitting gas runs diagonally from the northeast to south central in 
the map. A less bright irregular structure extends to the northwest from the south. 

The p-v diagrams (Figure [30]) show that the gas observable in the [Ne II] line has a steep 
velocity gradient. The [Ne II] emission near B2 is blueshifted by a few km s~^ with respect to 
the neutral cloud velocity. Over much of Bl, the bulk of the neon emission arises in redshifted 
material with the westernmost gas ranging from 0-20 km s~^ to the red of the molecular line 
center. The morphological differences between the radio and [Ne II] maps make it difficult to 



associ ate the kinematic and morphological structures. The RRL observations of iGaray et al. 



( I1994J ). which are not affected by extinction, show that the bulk of the ionized material is 
redshifted with respect to the molecular cloud. 



4.12. K3-50A (G70.3 +1.6) 



K3-50 is a four-com ponent (A-D. IWynn-Williamslll969l ) HII region complex at a distance 
of 8.7 kp c (|Harrislll975[) . The kinematic distance was computed using the galactic rotational 
curve of ISchmidtl (Il965l ) with the assumed Sun-to-GC distance of 10 kpc and a 250 km s~^ 
circular velocity at the position of the Sun (JRubin fc Turnerll 19691 : lRubinlll965l ). More recent 
estim ates of the distance to th e Galactic Center (szm7.94 kpc) imply that this distance is too 
large (JEisenhauer et al.l 120031 ). The UCHII regio n K3-50A appea r s to b e the youngest and 
dominates the emission at infrared wavelengths. iHofmann et al.l (120041 ) presented K' band 
bispectrum speckle interferometric data of K3-50A and their observations suggested that K3- 
50A is excited by a small cluster of massive to intermediate-mass stars. Radio continuum 
and hydrogen recombination line observat ions found that the ionized gas in K3-50A was 
undergoing a high- velocity bipolar outflow (Ide Free et al.lll994l ). The left panel of Figure 
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shows the radio continuum image (Ide Pree et al.lll994l ). There is extended emission both 
north and south of the dashed box that shows the extent of the [Ne II] map. 

The [Ne II] line emission in K3-50A peaks at the radio peak and has a remarkably 
similar distribution to that of the radio continuum within t he dashed box (Figure [31]) . Thi s 
is consistent with its relatively low extinction {Ak = 1.6, iMartin-Hernandez et al.ll2002l ). 
Most of the line emission comes from the core of the region. The integrated line map shows 
that the core of K3-50A consists of two sources, source 1 (SI) at the hne emission peak and 
source 2 (S2) ~ 1.5" southeast of SI. Line emission extends north and south from the central 
region. In the 12.8 /im continuum (Figure [3T|) . a more compact emission peak is slightly 
offset (~1") to the north of SI. 

The [Ne II] emission has a very broad velocity range, from -67 km s^^ to -8 km s~^ 
(Figure [5^ . SI, at (0,0), is present in all velocity channels. S2 is bright in the channels at 
the red end of the spectrum. A third component (S3) ~ 2" to the south of SI is present 
in the channels at central and blue-shifted velocities. All three cor nponents lie significantly 
to th e blue of the molecular cloud which has Vlsr = -24 km s~^ ( IShepherd fc Churchwell 
19961 ). S3 has the biggest velocity shift with a line center at ~-39 km s~^, with SI and S2 at 
~-35 km s^^. Position- velocity diagrams of K3-50A show the velocity differences between 
different sources (Figure [33]) . The position- velocity diagram for a north-south cut (top panel 
of Figure [33]) shows a velocity gradient while the region right around SI shows evidence for 
a barely resolved bipolar structure with a total velocity extent of 60 km s~^ (middle panel of 
Figure [55]) . Because the region is very extended north-south, line emission from the region 
fills our slit from end to end and emission from the northern and southern lobes are very 
faint. The emission indicated by the last four contours in Figure [31] is at least 100 times 
fainter than emission at the peak. Our map did not cover enough area to show the two more 
extended emission lobes. 



4.13. S106 

S106 (G76.4 -0.6) is a well-known bipolar HII region consisting of two lobes with the 
symmetry axis oriented at a position angle of ~30°. An equatorial gap is seen in all wave- 
lengths from optical to radio. Earlier investigators suggest that a dense, circumstellar disk 
surrounds a star and bl ocks the ionizing photons in the disk plane forming a biconical nebula 
( JBally fc Scovillelll982l ). A point source, pr esumably the excitation source of th e nebula, is 



l ocate d at the center of the equatorial gap (jPipher et al.|[l976l : [Felli et al.[[ 19841 ). [Felli et al. 



( jl984[ ) noticed that the bright radio continuum peaks in S106 tended to be distributed at the 
edges of the radio lobes and inferred that the inner volume was filled with gas with lower den- 
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sity. A radio continuum map from their observations is shown in Figure [SH Dense molecular 
clumps were identified on both sides of the central source of the HII region (ISchneider et al. 



2OO2I ) but there is no evidence of a smooth molecular disk. ISmith et al.l (1200 ll ) imaged the 
region at 3 to 20 /im. Some of their images (11.7 /im) show two dark lanes. One crosses 
the center of the region from the east to the west and divides the region into the north and 
south lobes. The other originates east of the region, crosses the south lobe and forms a small 
angle (~30°) with the first dark lane. 

The integrated [Ne II] line map and the radio continuum map are stunningly similar 
(Figure [M|) . The bipolar morphology and individual emission concentrations are clearly seen 
in both maps. The south lobe is characterized by a few bright clumps along the edge, while 
the north lobe contains more diffuse emission. The axis of the bipolar structure, which is 
indicated by the sharp edges of the south lobe, runs NNE-SSW. The main difference between 
the radio and [Ne II] maps is that the second dark lane seen in the infrared continuum is 
also seen in [Ne II]. 

The velocity structure shown in the p-v diagrams (Figure 1551) and the channel maps 
(Figure [3H]) is consistent with a picture of S106 as a nearly edge-on bipolar nebula. Both 
figures indicate that the south lobe is slightly blue shifted (~3 km s~^) with respect to the 
molecular cloud velocity while the material in the northern lobe is on-average red-shifted 
by ~4 km s~^. At the same time, individual intensity peaks show moderately narrow (~10 
km s~^ lines. Projection effects may lead to velocity shifts from one emission peak to its 
neighbors. 



4.14. NGC7538A 



NGC7538A belo ngs to an HII region-molecular cloud complex at a distance of ~3.5 kpc 
(JHanson et al.ll2002l ). In t he 6 cm radio continuum im age, it is a compact, spherical HII 
region with a radius ~ 5" (jWood fc Churchwelll Il989bl ) (Figure ETl left), corresponding to 
~19,000 AU. Molecular line ob servations show th at molecular gas towar d NGC7538 has 
a LSR velocity —56.9 km s~^ (JDickel et al.l Il98ll ). iBloomer et al.l (119981 ) suggested that 
NGC7538A is formed by a stellar wind bow shock with a star moving from the northwest to 
the southeast away from the Earth into the molecular cloud. 



Since dust extinction toward the region is not huge {A^ ~ 1-3, I Werner et al .11 19 79 



Herter et al.lll98ll : ICampbell fc ThompsonI 119841 ) . our [Ne II] line map of NGC7538A (Fig- 
ure [371 right) resembles the radio continuum map closely except that the radio component 
which should be at the southwestern corner of the source A, ~(-5,-2), is not present in our 
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[Ne II] map. This is probably due to our poor spatial resolution. [Ne II] line channel maps, 
spanning a velocity range from -82 to -51 km s~^, are shown in Figure [381 A cavity can 
be seen in the most blue-shifted channels and in the integrated line map. The shape of 
the region changes from ring-like at the blue-shifted channels to core-like at the red-shifted 
channels. The channel maps in Figure [38] show that most of the ionized gas in NGC7538A 
is blue-shifted with respect to the molecular cloud. The limb-brightened appearance of the 
source in the integrated line strength and radio continuum maps (Figure [37D might lead to 
the conclusion that the source is a limb-brightened spherical shell. However, the pattern of a 
small central source at the molecular cloud velocity and a larger, shell-like structure at bluer 
velocities, combined with the "A" and "7" structures seen in the pv cuts through the source 
(Figure [39l) argue against this picture. The kinematics and morphology are more consistent 
with a tangential flow along a shell with the symmetry axis close to the line of sight and 
with the vertex pointing away from the observer. 



4.15. W3A+B 



W3 belongs to an ex tensive HII region-mol ecular cloud complex located in the Perseus 
spiral arm (D = 2.0 kpc. iHachisuka et al.ll2006l ). Aperture synthesis observations identified 
four compact /ultracom pact radio continuum sources in a 3 arc minute diameter region (A-D, 
Wvnn-Williamsl 1971). Infrared ina aging at 1.65 to 20^m revealed nine objects (IRS 1-9) 
in W3 (IWynn- Williams et al.lll972l ). four of which are associated with the compact radio 
continuum sources: IRSl and IRS2 with the circular source W3A; IRS3 with W3B and 
IRS4 with W3C. Earlier authors have used the "broken-shell" morphology and the small 
velocity gradients seen in hyd rogen recombination lines toward W3A t o argue for a picture 



of this region as a blister flow (JDickel et al.lll983l: 



characterized as an emerging blister flow (JTieftrunk et al 



Tieftrun 



19971 ) 



k et al.lll997r) . W3B has also been 



W3A has the morphology of a limb-brightened broken shell in our [Ne II] map (Fig- 
ure [ID]), in close agreement with the distribution of the radio continuum emission. Most of 
the line emission forms a crescent open to the south with a horizontal structure crossing the 
center of the incomplete shell. An irregular emission cavity is present inside the crescent. 
Two finger-like structures extend from the southeast corner of the ring to the south. If we 
ignore these finger-like structures and another small lobe of emission on the northwest edge 
of the shell, the region closely resembles a half-ring. Such an image becomes clear in the 
channel maps of the region (Figure [H]). The half-ring morphology is present in channels 
over a broad velocity range. The sharp outer edge of the ring is evident. The emission level 
drops rapidly to the background level within a few arcseconds. 
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The shell-like morphology of W3A might lead one to expect the source to have the 
kinematics of an expanding shell. However, p-v diagrams (Figure H2l) of the region do not 
support this picture. We do not see the two components toward the center of the source 
which, in an expanding shell, would arise from the front and back sides of the shell and 
would gradually merge into one component at an intermediate velocity at the edge of the 
shell. Instead, we see a very broad line (>25 km s~^) along the bright arc-like rim of the 
region becoming narrow blue-shifted and red-shifted lines toward the center of the region. 
The change of the line width and the center velocities of the lines are abrupt. A surface 
flow, rather than an expanding shell, can better explain the observed p-v structure. The fact 
that both p-v cuts show the same pattern of a broad feature at the bright limb followed by 
a double line within the cavity implies that the vertex of the flow lies between the two cuts. 
An examination of the -53.5 km s~^ and -26.8 km s~^ channel maps in Figure HH supports 
this picture. The emission in the southern part of W3A comes from a narrower velocity 
range and therefore is not part of the surface flow seen in the north. 

The compact HII regions W3B and W3A have similar central velocities and line widths. 
The overall [Ne II] line emission distribution in W3B (see left panels in Figure H3ll looks 
like an irregular shell. The line proflle at most positions in the region can be flt by a single 
Gaussian. W3B contains a low surface brightness cavity, which can be seen in both the 
integrated line map and the channel maps (Figure l44l) of the region. The total shift in 
velocity centroid in W3B is relatively modest, from ~-42 km s^^ to -34 km s~^ and there is 
an apparently random element to the velocity variations at the 2-3 km s~^ level. Nonetheless, 
both the channel maps (Figure IH]) and the p-v cuts (Figure B31) are consistent with tangential 
flow along the ionized walls of a cavity. The channel maps show emission from the central 
region primarily at the least negative velocities, while the shell is most prominent in the 
bluer channels. The p-v diagrams also show the reddest velocities at positions crossing the 
center of the cavity and the bluest velocities along the edge. 



5. Discussion 

We have presented high spatial and spectral resolution mapping of the 12.8 ^m [Ne II] 
line toward a sample of 18 compact and ultracompact HII regions, where we count distinct 
HII regions within clusters as separate sources. In our disc ussion, we make use of previously 
pubhshed results for two of the 18 regions: G29.96-0.02 Jzhu et all boo si) and W51 IRS 2 



(JLacy et al.l 120071 ). We also include in the discussion one source (Mon R2. |Jaffe et al.ll2003l ) 



for which no results are presented in the current paper, bringing our sample to 19 sources. 
In order to discuss possible interpretations of the observed source morphologies and 
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kinematics, we first sketch p-v diagrams of several types of models that have been suggested 
for UCHII regions. We then compare the observed morphologies and kinematics to these 
sketches. We emphasize that the sketches are only meant to be qualitative descriptions of 
generic models. They are not derived from hydrodynamic simulations. 



5.1. Qualitative Models 

To facilitate the discussion of the observed gas motions, we first consider qualitative 
kinematic models of UCHII regions. Sketches of several types of models of UCHII regions 
are shown in Figures | 



A. (Figure H5|) An uniformly expanding shell source appears as a limb-brightened disk. 
Its p-v diagram is ring-like, centered on the molecular cloud velocity, Yamb- 



B. (Figure HS]) An expanding filled sphere appears as a limb-darkened disk. Its p-v 
diagram is a filled ellipse centered at Yamh- 



C. (Figure 1461) The appearance of a bow-shock source, in which a wind from a star 
moving supersonically through a molecular cloud sweeps up a paraboloidal shell, with the 
motion of the gas in the s hell resulting from the combined momenta of the stellar wind 



and the swept-up material (|Mac Low et al.lll99ll : Ivan Buren fc Mac Lowlll992l : IWilkinI 12000 



Zhu et al.l l2005l : [Arthur fc Hoard 120061 ) depends on the viewing angle. If viewed from the 
side {9 = 90°, left panel), two velocities are seen at each offset, on opposite sides of Vamb- 
The line is broad in the head, and the two velocities approach each other in the tail. If 
viewed from the tail {6 = 180°, right panel), the source would not be obviously cometary in 
appearance, and the velocity would vary from V* at the center to Vamb at the edges. From an 
intermediate angle [6 = 135°, middle panel), we see a p-v pattern like those seen in G29.96, 
but offset in velocity. A broad line, centered between Vamb and V^ is seen at the head of the 
comet. Of the two velocities seen at 6* = 90°, the V < Vamb branch is foreshortened and may 
blend into the head (see Figure H]). The V > Vamb branch rises toward V^ then falls back 
toward V = Vamb- The line centroid is at V > Vamb- The p-v diagrams would be mirrored 
about Vamb if the star were moving toward the observer {9 < 90°). 

D. (Figure HTj) A paraboloidal shell can also occur in which the motion of the gas in 
the shell is predominantly due to the pressure gradient along the shell, wh ich causes an 
acceleration of the gas toward the tail (IZhu et al.ll2005l : [Arthur &: Hoard 120061 ). In this case, 
the paraboloidal shape of the shell could result either from a relatively slow motion of the 
star through the molecular cloud or from a density gradient in the cloud. The resulting 
appearance and p-v diagrams are similar to those of the bow-shock model, but the p-v 
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diagrams are offset toward negative velocities, with the velocity centroid at V < Vamb (for 
9 > 90°). This is the offset seen in G29.96, indicating that it belongs to this category. 



E. (Figure H51) A filled blister flow could result if the wind from the ionizing star is too 
weak to preve nt the ionized gas from fiow'ing away from the surface of the molecular cloud , 
past the star (JYorke et al.lll983l : IComeronI 119971 : iHenney et al.ll2005l : [Arthur fc Hoard l2006l ). 
Such a region is formed when a density gradient is present in the natal molecular cloud. In 
this case, the source morphology and p-v diagrams would depend on the shape of the blister 
cavity in the molecular cloud, but blister flows are largely axial and a broad line would be 
likely to be seen at all positions. For 9 > 90°, the line centroid is at V < Vamb- The largest 
velocities would be seen toward the tail, or at the center for 9 ~ 180°. 

A number of sources in our sample either have confused kinematics or have kinematics 
for which we have no simple qualitative models: Several sources show ring, horseshoe, or 
bipolar morphologies, without the kinematic signatures of shells or paraboloidal flows. Sev- 
eral regions appear to be clusters of overlapping sources, making the interpretation of their 
kinematics difficult. Several sources have distinctly different morphologies in [Ne II] and the 
radio continuum, suggesting that variable extinction is distorting their appearance in the 
infrared. These sources also have confusing kinematics. And several sources are too small 
for their [Ne II] kinematics to have been resolved. 



5.2. Observed Source Morphologies 

The extraordinarily good match between the radio and [Ne II] morphologies in the ma- 
jority (12 of 16) of the resolved sources indicates that most neon is in Ne"*", that the gas 
density is almost everywhere below the critical density of the 12.8 fim line, and that the 
centimeter-wave radio continuum emission is for the most part optically thin. Extinction 
is the most likely cause of the strong differences between radio continuum and [Ne II] mor- 
phologies in the other four sources. Some of the differences could also be due to incomplete 
UV coverage in the radio snapshot observations. And, differences in sensitivity between 
[Nell] and radio continuum observations may also play a role. 

In many cases, radio interferometer observations show that UCHII regions are accompa- 
nied by substantial radio continuum flux from extended, low-surface brightness ionized gas 
(JKurtz et al.lll999l : iKim fc Kodl200ll ). These results imply that UCHII regions may exist in 
asymmetric regions where the ionizing radiation reaches dense gas in one direction and more 
extended low density gas in the other direction. Our scanning technique does not permit us 
to find these extended haloes in [Ne II] due to limited sky area we can cover and the sky 
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background subtraction method we use, but the frequent presence of emission halos may 
have imphcations for the interpretation of the source kinematics. 



5.3. Observed Source Kinematics 

Of the 12 sources that are resolved and not badly confused by extinction, nine show 
clear evidence of bulk motions, which are usually not symmetric about the molecular cloud 
velocities. These nine sources all show strong evidence of tangential flows along the surfaces 
of shells. Three of these sources, G29.96 (Figures [3] and HI), G33.92 (Figure [13]), and 043.89 
(Figure [T71) . show, with varying degrees of clarity, the "A" and "7" patterns in p-v diagrams 
that are similar to the patterns established in the surface flow models (types C and D). 
In these three cases, the patterns indicate that the terminal velocities are blueshifted with 
respect to the molecular cloud velocities. Along the symmetry axis, the velocities cross the 
molecular cloud velocity from blue to red and back to blue going from the vertex into the 
source for 029.96 and 033.92. In 043.89, the velocity approaches the cloud velocity from the 
blue side but does not cross it before moving back to the blue. The kinematic behavior we 
observe in all three sources is very similar to what we see in the models of 029.96 (Figures [3] 
and HI) , which are tilted with the vertex away from the observer and have an inclination 
of the symmetry axis of about 45° with respect to the line of sight. With their vertices 
pointed away from the observer, the blueshift relative to Vamb means that these sources are 
of kinematic type D (pressure-driven surface flows; Figure Wl\ . 

The three sources in our sample whose [Ne II] and radio continuum morphologies most 
closely resemble limb-brightened spherical shells do not have kinematic patterns consistent 
with a closed-shell geometry. The kinematic pattern s for three of the sources, NOG 7538A 



(Figure [391). W3B (Figure [l3l). and W51 IRS2, fsee [Lacv et all [20071 ) . are consistent with 
motion of the ionized gas along roughly parabolic shells with vertices pointing almost directly 
away from the observer. In the position-velocity plots, the most positive velocities, usually 
close to the velocity of the dense molecular gas, are toward the vertex while the velocities 
along the bright rim are shifted to the blue. The presence of this pattern is largely inde- 
pendent of the position angle of the position-velocity cut through the center of the source. 
These sources are of kinematic type D with 9 ~ 180°. 

The remaining sources with understandable kinematics also have a strongly limb-brightened 
appearance. The kinematics revealed by the [Ne II] sp ectroscopy show that 030.54 (Fig- 



ure[II]), W3A (Figure[l2]), and Mon R2 ([jaffe et al.ll2003l ) are also likely to be limb-brightened 



parabolas with significant flow of the ionized gas along the ionized shell. For these three 
sources, however, the kinematics indicate that the symmetry axes lie near the plane of the 



-22- 



sky. The sources show broad hnes at the source vertex and then double hues farther along 
the symmetry axis with the two components appearing at the extreme velocities of the line 
seen at the vertex. Along the edges of each source, the lines are single and significantly 
narrower. These sources are of either kinematic type D or type C, with 9 ^ 90°. 

Four sources have very similar [Ne II] and radio morphologies, but do not fit into any of 
the kinematic categories we have discussed. S106 (Figure [Ml) has a clear bipolar morphology, 
but confused kinematics. G45.12N (Figure fITh has a cometary morphology, but does not 
have obvious cometary kinematics, of either type C or type D. G45.12SW (Figure 1221) has 
a ring or shell morphology, but does not have the kinematic signature of a shell. K3-50A 
(Figure [31]) appears to be a cluster of sources, which are too blended to allow us to sort out 
their kinematics. The peaks in the IRTF map of [Ne II] in W51 IRS2 (Figure [2S]) are also 



blended. The Gemini map of [Lacy et al.[ (|2007[ ) separates these sources but does not resolve 



th eir kinemat i cs. H owever, the kinematics of the extended ionized gas in W51 IRS2 shown 



m 



Lacy et al.[ (I20071 ) look much like a type D surface flow with 9 ~ 180°. 



Complex morphologies and kinematics in the four sources with substantially different 
[Ne II] and radio continuum morphologies prevent a coherent explanation for gas motion in 
these sources. However, it is clear that systematic gas motions are present. In Gil. 94, there 
is a component of the [Ne II] emission at the same velocity as the ambient cloud (Figure [8]), 
but there is also a substantial amount of emission at many positions in a red wing or in 
an identifiable red component. In G45.45 (Figure [25]) . the kinematic pattern is confusing, 
and the most prominent feature in [Ne II] is not apparent in the radio continuum. G61.48 
(Figure [29]) is the source where the [Ne II] morphology least resembles the distribution of the 
radio continuum emission. The [Ne II] emission we see may therefore not be sampling the 



bulk of the gas, but the dominance of red-shifted gas in the RRL observations (JGaray et al. 



1994[ ) is consistent with the conclusion from the [Ne II] -radio continuum morphology dif- 
ference, namely that the source is on the back side of the molecular cloud and the ionized 
material is fiowing away from the neutral core. 



5.4. Comparison of Observed Sources with Qualitative Models 

The kinematic types that we have assigned to observed sources are given in Table [2] 
Examination of this list shows that none of the observed sources fits into category A or B, 
the classical expanding shell and sphere models of HII regions. None has a ring-like p-v 
diagram of an expanding shell, and the sources with broad lines are either too small for their 
morphologies to be apparent or have confused morphologies. In addition, none of the sources 
fits into category E, the classical filled blister fiow with the ionized gas fiowing away from 
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the ionization front. 

Even more remarkably, in spite of the frequency of bow-shock morphologies, none of 
the observed sources clearly fits into category C, a bow shock moving with a star through a 
molecular cloud. Because they are nearly edge-on, Mon R2, G30.54, and W3A could fit into 
either category C or category D. The other sources with cometary kinematic signatures all 
have p-v patterns characteristic of sources tipped away from the observer {9 > 90°) with their 
velocity centroids blueshifted relative to the neighboring molecular material, the signature 
of a pressure-driven surface flow. 

The preponderance of sources with their tails tipped toward the observer {9 > 90°) 
probably results from a combination of extinction and selection effects. The extincti on at 
12.8 /im is relatively low fAi>> «/Av^0.05. lRieke fc LebofskvlEgSsI : [Rosenthal et allboooh . but 
the column density through the kind of massive core tha t harbors massive yo ung stars and 
UCHII regions can be substantial (Ay ^ a few hundred, iBeuther et al.l 120001). The density 



in such cores drops rather steeply, with power-law exponents of 1.6-1.8 (JMueller et al.l 12002 



Beuther et al.l 120021 ). If UCHII regions are not in the centers of cores, they should therefore 
be visible in [Ne II] when the exciting stars lie between the front side and 9 ^ 90°. Our 
survey was somewhat biased in favor of sources lying on the front sides of molecular clouds 
by including sources known from optical and near-infrared observations (NGC 7538, W3, 
and Mon R2). In addition, we were forced to abandon efforts to observe several UCHII 
regions (G8.67-0.36, GIO. 10+0.74, GlO.62-0.38, G12.208-0.10, G33.50+0.20, G42.42-0.27, 
G45.47+0.05, G50.23+0.33) because they were undetectable or too faint to map. These 
selection effects are likely to explain the lack of sources with kinematics indicating that they 
lie on the back sides of their clouds. 

More than half of the observed sources fit into category D, a flow along the surface of a 
paraboloidal shell, which accelerates relative to molecular cloud as it moves away from the 
apex of the paraboloid. Although we have not made a hydrodynamic model of such a flow, 
we know of no driving mechan ism for the acceleration other than the pressure gradient in the 
ionized gas. IZhu et al.l (120051 ) found that the pressure gradient acceleration would be small 
if the ionized and swept-up molecular gas were coupled so both had to be accelerated. But 
we have estimated the acceleration if the ionized gas is allowed to slip past the molecular 
shell, and we find tha t speeds ~20-30 km s~^ can occur. The hydrodynamic models of 



Arthur fc Hoard (120061 ) include the essential effects present in the combination of blister and 
bow shock models: de n sity g radients in the neutral medium, stellar winds, and stellar motion. 
The [Arthur fc Hoard (120061 ) models do an excellent job of matching the range of source 
morphologies and reproduce the most significant kinematic signatures using reasonable values 
for density gradients, mass loss rates, and stellar velocities. One critical area where they 



-24- 



fall short is in explaining the width of the hnes at t he source vertex . The observed effect 
is most evident in W3A (Figure H2l) and in Mon R2 (jjaffe et aLll2003l ) where the symmetry 
axes are near the plane of the sky. The splitting of the broad line into two narrow lines at 
the velocity extremes of the broad line as one crosses the shell to the inside implies that the 
gas accelerates to ~ ±10 km s~^ within a small solid angle (less than a steradian, as seen 
from the point of view of the ionizing star) at the vertex. This, in turn, implies a significant 
pressure gradient. 

The [Ne II] kinematics of our sample unite a large and diverse group of UCHII and 
compact HII regions into a single empirical picture: The HII regions coexist with dens e 
and massive molecular cores (JBeuther et al.ll2002l : IShirley et al.ll2003l : ISridharan et al.ll2002l ). 
The ionized regions have the form of partial, roughly parabolic shells with their open sides 
usually pointing away from the centers of the dense neutral condensations and gas flows 
tangentially along these shells. 

The evidence that very many UCHII regions are dense shells with tangential flows 
away from the parent molecular cores has implications for the survivability of the molecular 
cores. Earlier examinations of the erosion of dense cores by external O stars did not take 
the compression of neutral material by stellar wind shocks into account. These efforts led to 
estimates that when as little as 4% of the core material forr ned into stars with a Salpeter IMF, 
these stars could erode the remaining core entirely away (jWhitworthl Il979l ) . It has already 
been pointed out, however that the dense shells formed by stellar wind pressure would g reatly 
inhibit the process of erosion (j Chur chwelll 1 1 9 9 ll : iLumsden fc Hoard Il999l : lHoardl2006l ). Our 
results indicate that most if not all OB stars breaking out of dense cores form such dense 
shells at their boundaries. 
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Fig. 1. — 2 cm continuum map (left, iFey et al.l Il995l . restoring beam: 0."56 x 0."49) and 
integrated [Ne II] line map (right) of G29.96 -0.02. The (0,0) position marks the peak 
emission in the [Ne II] map. The line map has been cross-correlated with the continuum map 
and shifted to match. Contours are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% 
of the peak value in each map. The peak [Ne II] surface brightness is 1.0 ergs cm~^ s~^ sr~^. 
The total [Ne II] flux is 6.3x10"^° ergs cm'^ s'^ 
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Fig. 2. — Channel maps of [Ne II] line observations for G29.96 -0.02. The individual 
maps sample a single velocity slice in the data cube and therefore represent a width equal 
to the spectral resolution, or about 4 km s~^. Contours are drawn at 70%, 50%, 35%, 
25%, 17.5%, 12.5%, 9% and 6% of the peak value of all channels. The peak value is 
50.0 ergs cm~^ s^^ sr^^ (cm^^)^^. The molecular cloud velocity is Vlsr = 98 km s~^. 
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Fig. 3. — Position- Veloc ity (P-V) diaKra ms of G29.96 -0.02 [Ne II] observations (column 2) 
and a bow shock model (jZhu et al.ll2005l . column 4). The straight lines through the observed 
(column 1) and model (column 3) integrated line flux maps show the position of each cut. 
Dashed lines in p-v diagrams show the observed ambient molecular material velocity. Dotted 
lines in the model p-v diagrams show the ambient material velocity in the model. The contour 
maps are drawn as in Fig{T]with the peak value ~1.0 ergs cm~^ s^^ sr~^ for the line maps 
of observations. 
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Fig. 4. — - Positi on- velocity diagrams of G29.96 -0.02 [Ne II] observations (left) and the 
Zhu et al.l (120051 ) bow shock model (right). The straight line in the integrated line map 
above each position- velocity diagram shows the location of the cut. The dashed lines in the 
p-v diagrams show the ambient molecular material velocity. Dotted lines in the model p-v 
diagrams show the ambient material velocity in the model. The contours are drawn as in 
Fig{T]with the peak value ~1.0 ergs cm^^ s^^ sr^^ for the line maps of observations. 
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Fig. 5. — 6 cm continuum map (left, IWood fc Churchwelll Il989bl . restoring beam: 0."38 x 
0."38) and integrated [Ne II] line map (right) of G5.89 -0.39. The [Ne II] line map has 
been cross-correlated with the radio map and the (0,0) position has been shifted to the best 
match. The (0,0) position is the location of the emission peak in the line map. Contours are 
drawn at 95%, 90%, 80%, 70%, 60%, 50%, 40% ,30%, 20% and 10% of peak value for the 
continuum map and at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value 
for the line map. The peak [Ne II] surface brightness is 1.08 ergs cm^^ s^^ sr^^. The total 
[Ne II] flux is 2.26x10-1° ergs cm'^ g-^. 
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Fig. 6. — Channel maps of G5.89 -0.39 [Ne II] line observations. Contours are drawn at 70%, 
50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value of all channels. The peak value 
is 9.15 ergs cm"^ s~^ sr~^ (cm^^)^^. The molecular cloud velocity is Vlsr = 9 km s~^. 
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Fig. 7. — 2 cm continuum map (left, IWood fc Churchwelj Il989bl . restoring beam: 0."41 x 
0."41) and integrated [Ne II] line flux map (right) of Gil. 94 -0.62. The (0,0) is the location 
of the peak [Ne II] line emission. The radio continuum peak, mark ed with an "X" , i s offse t 
1.17" west and 1.67" north from the line emission peak according to lDe Buizer et al.l (120031 ). 
Contours are drawn at 95%, 90%, 80%, 70%, 60%, 50%, 40% ,30%, 20%,15%, 10% and 5% of 
the peak flux density for the continuum map and at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% 
and 6% of the peak value for the line map. Five separate components (A, B, C, D and E) are 
indicated in the line map. The peak [Ne II] surface brightness is 0.073 ergs cm~^ s^^ sr^^. 
The total [Ne II] flux is 6.3 xlO"" ergs cm'^ s"^ 
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Fig. 8. — [Ne II] line profiles across Gil. 94 -0.62 (A, B and C). Spectra are averages over 1.5" 
xl.5". Central coordinates of the averaging areas are shown in parentheses. The brightness 
has units of erg cm^^ s^^ sr^^ (cni^^)^^. The pixel size along the spectral direction is 
0.0025 cm^^ (~0.95 km s^^). The molecular cloud velocity is Vlsr = 39 km s^^. 
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Fig. 9. — 2 cm continuum map (left, IWood &: Churchwelj Il989bl . restoring beam: 0."41 x 
0."41) and integrated [Ne II] line map (right) of G30.54 +0.02. The [Ne II] line map is cross- 
correlated with the continuum map and shifted to the best match. The (0,0) position marks 
the peak emission in the line map. Contours are drawn at 95%, 90%, 80%, 70%, 60%, 50%, 
40% ,30%, 20% and 10% of the peak value for the continuum map and at 70%, 50%, 35%, 
25%, 17.5%, 12.5%, 9% and 6% of the peak value for the line map The peak [Ne II] surface 
brightness is 0.086 ergs cm" 
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The total [Ne II] flux is 4.8x10"^^ ergs cm'^ s'^ 
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Fig. 10. — [Ne II] line profiles for G30.54 +0.02. Spectra are averages over 1.5" x 1.5". Central 
coordinates of the averaging areas are shown in parentheses. The brightness has units of 
erg cm~^ s~^ sr~^ (cm~^)~^. The pixel size along the spectral direction is 0.0025 cm~^. The 
molecular cloud velocity is Vlsr = 48 km s^^. 
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Fig. 11. — Position-velocity diagrams of G30.54 +0.02. Contours in line flux maps 
are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value 
(0.086 ergs cm~^ s~^ sr^^). Dashed lines in p-v diagrams show the ambient molecular 
material velocity. 



-42- 



00' 52 05" 



g5l'55"h 




RIGHT ASCENSIOW 




-5 



Fig. 12. — 6 cm continuum map (left. lFey et al.lll992l . restoring beam: l."6 x l."4, positions 
in B1950 coordinates) and integrated [Ne II] line flux map (right) of G33.92 +0.11. The box 
drawn with dashed lines in the continuum map indicates the approximate area covered by 
the [Ne II] line observations. The (0,0) position in the line map is the location of the peak 
[Ne II] line emission. Contours were drawn at 100%, 90%, 80%, 70%, 60%, 50%, 40%, 30%, 
20%, 10% and 5% of the peak value for the 6 cm continuum map and at 70%, 50%, 35%, 25%, 
17.5%, 12.5%, 9% and 6% of the peak value for the [Ne II] line map. The peak [Ne II] surface 



brightness is 0.24 ergs cm^ 
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The total [Ne II] flux is 2.7x10"^° ergs cm^^ g-^. 
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Fig. 13. — Position- velocity diagrams of [Ne II] observations toward G33.92 +0.11. Contours 
in line flux maps are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak 
values (0.24 ergs cm~^ s^^ sr^^) of the maps. Dashed lines in p-v diagrams show the ambient 
molecular material velocity. 
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Fig. 14. — [Ne II] line profiles for G33.92 +0.11. Spectra are made by averaging spectra 
over 1.8" x 1.8". Central coordinates of the averaging areas are shown in parentheses. The 
brightness has units of erg cm^^ s^^ sr^^ (cm^^)^^. The pixel size along the spectral direction 
is 0.0025 cm-^ 
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Fig. 15. — Channel maps of [Ne II] line observations for G33.92 +0.11. Contours 
are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value 
(5.9 ergs cm^^ s^^ sr^^ (cm^^)^^). The molecular cloud velocity is Vlsr = 108 km s~^. 
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Fig. 16. — 2 cm continuum contour map (left. IWood fc Churchwel]||l989bl . restoring beam: 
0."41 X 0."41) and integrated [Ne II] line flux map (right) of G43.89 -0.78. The line map 
has been cross-correlated with the continuum map and shifted to match. The (0,0) position 
is the location of peak emission in the line map. Contours are drawn at 95%, 90%, 80%, 
70%, 60%, 50%, 40%, 30%, 20%, 15%, 10% and 5% of the peak value for the radio map 
and at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value for the line map. 
The peak [Ne II] surface brightness is 0.22 ergs cm~^ s~^ sr~^. The total [Ne II] flux is 
7.3x10^^^ ergs cni^^ s^^. 
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Fig. 17. — Position- velocity diagrams of [Ne II] line observations (two columns on the left) 
and a bow shock model (two columns on the right) of G43.89 -0.78. Contours in line flux 
maps are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak values 
(0.22 ergs cm~^ s^^ sr^^) in the maps. See text for a description of the model. The dashed 
lines in the p-v diagrams show the ambient molecular cloud velocity. Dotted lines in the 
model p-v diagrams show the ambient material velocity in the model. 
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Fig. 18. — Dust 12.8 /iin continuum map (left) and integrated [Ne II] line flux map (right) of 
G45.07 +0.13. The fainter northern source in the left figure is a hot molecular core candidate. 
The (0,0) position marks the peak [Ne II] line emission. Contours are drawn at 70%, 50%, 
35%, 25%, 17.5%, 12.5%, 9% and 6% of peak values in both maps. The peak continuum 



flux density is 0.23 ergs cm 
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^ The total [Ne II] flux is 7.8xl0~^^ ergs cm'^ s~^ 
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Fig. 19. — 2 cm continuum map (left. IWood &: Churchwel]||l989bl . restoring beam: 0."41 x 
0."41) and integrated [Ne II] line map (right) of G45.12 +0.13. The (0,0) position is the 
location of the [Ne II] line emission peak. The contours are drawn at 95%, 90%, 80%, 70%, 
60%, 50%, 40% ,30%, 20%,15%, 10% and 5% of peak value for the radio continuum map 
and at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9%, 6%, 4%, 2% and 1% of the peak value for 
the [Ne II] map. The peak [Ne II] surface brightness is 1.2 ergs cm~^ s~^ sr~^. The total 
[Ne II] flux is 4.5x10"-*^° ergs cm~^ s~^ 
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Fig. 20.— Channel maps for G45.12 +0.13. Contours are drawn at 70%, 50%, 35%, 25%, 
17.5%, 12.5%, 9%, 6%, 4%, 2% and 1% of the peak value of all channels. The peak value is 
11.1 ergs cm~^ s~^ sr~^ (cni~^)~^. The molecular cloud velocity is Vlsr = 59 km s~^. 
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Fig. 21. — Position- velocity diagrams of G45.12 +0.13N. Contours in line flux maps 
are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak values 
(1.23 ergs cm~^ s~^ sr~^) in the maps. Dashed lines in p-v diagrams show the ambient 
molecular material velocity. 
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Fig. 22. — Position- velocity diagrams of G45.12 +0.13SE. Contours in line flux maps 
are drawn at 70%, 50%, 35%, 25%, f7.5%, f2.5%, 9% and 6% of the peak values 
(0.10 ergs cm^^ s~^ sr~^) in the maps. Dashed lines in p-v diagrams show the ambient 
molecular material velocity. 
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Fig. 23. — 6 cm continuum map (left. IWood fc Churchwel]||l989bl . restoring beam: 0."38 x 
0."38) and integrated [Ne II] line map (right) of G45.45 +0.06. The [Ne II] line map has 
been cross-correlated with the radio map and the (0,0) position has been shifted to the best 
match. The (0,0) position is the location of the emission peak in the line map. Contours are 
drawn at 95%, 90%, 80%, 70%, 60%, 50%, 40% ,30%, 20% and 10% of peak value for the 
continuum map and at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value 
for the line map. The peak [Ne II] surface brightness is 0.16 ergs cm~^ s~^ sr^^. The total 
[Ne II] flux is 1.2x10"^° ergs cm'^ g-^. 
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Fig. 24. — Channel maps for G45.45 +0.06 [Ne II] line observations. Contours are drawn 
at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value of all channels. The 
peak value is 4.2 ergs cm~^ s~^ sr~^ (cm~^)~^. The molecular cloud velocity is Vlsr = 59 
km s~^. 
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Fig. 25. — 2 cm continuum map (left. lGaume fc Mutellll987l . resolution: l."2) and integrated 
[Ne II] line flux map (right) of W51 IRS2. Contours in the radio map were drawn at 90%, 
80%, 70%, 60%, 50%, 40%, 30% ,20%, 10% and -10% of the peak value. Contours in the 
[Ne II] map are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of peak line 
emission at (0,0). The peak [Ne II] surface brightness is 0.66 ergs cm~^ s~^ sr~^. The total 
[Ne II] flux is 4.2x10"^° ergs cm'^ s"^ 
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Fig. 26. — Channel maps of W51 IRS2 [Ne II] line observations. Contours are drawn at 70%, 
50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value of all channels. The peak value 
is 13.2 ergs cm^^ s^^ sr^^ (cm^^)^^. The molecular cloud velocity is Vlsr = 58 km s~^. 
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Fig. 27. — Position- velocity diagrams of W51 IRS2 [Ne II] line observations. Contours in 
line flux maps are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of peak values 
(0.66 ergs cm~^ s~^ sr~^) of the maps. Dashed lines in p-v diagrams show the ambient 
molecular material velocity. 
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Fig. 28. — 1.3 cm continuum map of G61.48 +0.09 (iGomez et al.l Il995l . resolution ~ 4"). 
Contours were drawn at 90%, 75%, 60%, 45%, 30% ,20%, 10%, 5% and -5% of the peak 
value. The box drawn with dashed lines indicates the approximate area covered by our 
[Ne II] observations. 
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Fig. 29. — Integrated [Ne II] line flux map of G61.48 +0.09. The (0,0) position is the location 
of the [Ne II] line emission peak. Contours are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 
9% and 6% of the peak value. The peak [Ne II] surface brightness is 0.052 ergs cm^^ s^^ sr^^. 
The total [Ne II] flux is 2.2x10"^° ergs cm^^ g-i_ 
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Fig. 30. — Position-velocity diagrams of G61.48 +0.09 [Ne II] line observations. Contours in 
line flux maps are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of peak values 
(0.052 ergs cm~^ s~^ sr~^) in the maps. The spatial resolution of the map is ~ 4". Dashed 
lines in p-v diagrams show the ambient molecular material velocity. 
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Fig. 31. — The 2.0 cm continuum map of K3-50A (left) (Ide Free et al.l Il994j . resolution~ 
l."3). The (0,0) position marks the peak continuum emission. First positive and negative 
contours were drawn at the 3cr flux level. Other contours were drawn at 1.4, 2, 2.8, 4, 8, 
16, 32, 64, 128 and 256 times the 3cr flux level. The box drawn with dashed lines indicates 
the approximate area covered by our [Ne II] observations. The 12.8 /im dust continuum 
map (top right) and integrated [Ne II] line flux map (bottom right). For both maps, the 
(0,0) position is the location of the [Ne II] emission peak. Three crosses show the approx- 
imate positions of SI, S2, and S3 which are [Ne II] emission concentrations identified in 
the channel maps (Figure [32]). Contours are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 
9%, 6%, 4%, 2% and 1% of the peak values in both maps. The peak dust continuum 
flux density is 0.82 ergs cm~^ s~^ sr~^ (cm^^)^^. The peak [Ne II] surface brightness is 
0.34 ergs cm^^ s^^ sr^^. The total [Ne II] flux is 1.0x10"^° ergs cm~^ s~^. 
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Fig. 32. — Channel maps of K3-50A [Ne II] line observations. Contours are drawn at 70%, 
50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value of all channels. The peak value 
is 5.1 ergs cm^^ s^^ sr^^ (cm^^)^^. The molecular cloud velocity is Vlsr = -24 km s~^. 
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Fig. 33. — Position- velocity diagrams of K3-50A [Ne II] line observations. Contours in line 
flux maps are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of peak values 
(0.34 ergs cm~^ s~^ sr~^) of the maps. Dashed lines in p-v diagrams show the ambient 
molecular material velocity. 
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Fig. 34. — 1.3 cm continuum map (left. iFelli et al.lll984l . resolution^ 2. "3) and integrated 
[Ne II] line flux map (right) of S106. Contours in the continuum map were drawn at each 
5% step from 5% to 95% of the peak value. The box drawn with dashed lines indicates 
the approximate area covered by the [Ne II] observations. Contours in the [Ne II] line map 
are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak line emission. 
The spatial resolution of the radio map is ~ 2. "3. The (0,0) position is the location of a 
weak continuum source seen in the line observations. The peak [Ne II] surface brightness is 
0.23 ergs cm~^ s~^ sr~^. The total [Ne II] flux is 1.1x10"^ ergs cm~^ s~^. 
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Fig. 35. — Position-velocity diagrams of S106 [Ne II] line observations. Contours in [Ne II] 
line flux maps are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of peak values 
(0.23 ergs cm~^ s~^ sr~^) in the maps. Dashed lines in p-v diagrams show the ambient 
molecular material velocity. 
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Fig. 36. — Channel maps of S106 [Ne II] line observations. Contours are drawn at 70%, 50%, 
35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value of all channels. The peak value is 
3.6 ergs cm"^ s~^ sr~^ (cm~^)^^. The molecular cloud velocity is Vlsr = -1 km s~^. 
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Fig. 37. — 6 cm continuum map (left. IWood fc Churchwel]|ll989bl . restoring beam: 0."38 x 
0."38) and integrated [Ne II] line flux map (right) of NGC7538A. Contours in the continuum 
map were drawn at 95%, 90%, 80%, 70%, 60%, 50%, 40% ,30%, 20%, 15%, 10%, 5%, -5%, 
— 10%, —15% and —20% of the peak value. The box drawn with dashed lines indicates 
approximately the area covered by the [Ne II] observations. Contours in the [Ne II] line 
map are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value in 
the map. The radio map has a synthesized beam with FWHM~ 0."38 x 0."38. The (0,0) 
position is the location of the [Ne II] line emission peak. The peak [Ne II] surface brightness 
is 0.16 ergs cm~^ s~^ sr~^. The total [Ne II] flux is 2.7x10^^*^ ergs cm~^ s~^. 
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Fig. 38. — Channel maps of N7538 A [Ne II] line observations. Contours are drawn at 70%, 
50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value of all channels. The peak value 
is 3.1 ergs cm~^ s~^ sr~^ (cm~^)~^. The molecular cloud velocity is Vlsr = -57 km s~^. 
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Fig. 39. — Position-velocity diagrams of NGC7538 A [Ne II] line observations. Contours in 
line flux maps are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak 
values (0.16 ergs cm^^ s^^ sr^^) of the maps. Dashed lines in p-v diagrams show the ambient 
molecular material velocity. 
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Fig. 40. — Grey-sca le map of the radio 6 cm free-free continuum of W3 core (left) 



(JTieftrunk et al.l 119971 ) . The observations have a FWHM beam with a size of 0."71 x 0."56. 
The flux in the region ranges from 0.1 mJy to 40 mJy beam~^. The box drawn with dashed 
lines indicates the approximate area covered by our [Ne II] observations. Integrated [Ne II] 
line flux map of W3 A and B region (right). The (0,0) position is the peak location of 
integrated emission map. Contours are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% 
and 6% of the peak value. The peak [Ne II] surface brightness is 0.26 ergs cm~^ s~^ sr~^. 
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Fig. 41. — Channel maps of W3 A [Ne II] line observations. Contours are drawn at 70%, 
50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value of all channels. The peak value 
is 2.8 ergs cm~^ s~^ sr~^ (cm~^)~^. The molecular cloud velocity is Vlsr = -40 km s~^. 
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Fig. 42. — Position- velocity diagrams of W3 A line observations. Contours in line flux 
maps are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of peak values 
(0.13 ergs cm~^ s~^ sr~^) of the maps. Dashed lines in p-v diagrams show the ambient 
molecular material velocity. 
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Fig. 43. — Position-velocity diagrams of W3 B line observations. Contours in line flux 
maps are drawn at 70%, 50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of peak values 
(0.17 ergs cm~^ s"^ sr~^) of the maps. 
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Fig. 44. — Channel maps of W3 B [Ne II] line observations. Contours are drawn at 70%, 
50%, 35%, 25%, 17.5%, 12.5%, 9% and 6% of the peak value of all channels. The peak value 
is 5.2 ergs cm^^ s^^ sr~^ (cm^^)^^. The molecular cloud velocity is Vlsr = -43 km s^^. 
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Fig. 45. — Sketches of position-velocity diagrams of an expanding shell-like HII region (left) 
and an expanding spherical HII region (right). The expansion speed is V^. The molecular 
cloud velocity is Vamt- The dashed lines show the locations of the p-v cuts. The solid 
arrows in the morphology diagrams indicate the moving directions of ionized gas relative 
to ambient molecular material and the dashed arrows show the moving directions of stellar 
wind material with respect to the central stars. 
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Fig. 46. — Sketches of position-velocity diagrams of a stellar wind bow shock HIT region 
at three different viewing angles. Note that the velocity vectors are shown in the frame of 
reference of the ambient molecular cloud; the shell is not expanding in the moving frame of 
reference of the star. The molecular cloud velocity is Vamb and the ionizing star is moving 
supersonically at a speed of V* relative to the molecular cloud. The ram pressures of incoming 
cloud material and stellar wind material confine the ionized material in a compressed shell 
which moves with the star in the cloud. The dashed lines show the locations of the p-v cuts. 
The p-v diagrams in the middle and right panels are drawn for stars moving away from the 
observer (6 > 90°). The p-v diagrams would be mirrored about Vamb for stellar motions 
toward the observer (6 < 90°). 
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Fig. 47. — Sketches of position- velocity diagrams of a stellar wind pressure-driven HIT region 
at three different viewing angles. The molecular cloud velocity is Vamb and the star with a 
stellar wind is stationary with the cloud. The ionized gas is compressed by the stellar wind 
into a shell and moves along the shell under the influences of the density gradient and the 
stellar wind pressure. Voo is the exiting speed of the ionized gas. The dashed lines show the 
locations of the p-v cuts. The p-v diagrams are similar to those in Fig. HHl but are offset 
toward negative Ylsb. for the head of the shell tipped away from the observer {9 > 90°). 
The p-v diagrams would be mirrored about Vamb for stellar motions toward the observer 
{9 < 90°). 
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Fig. 48. — Sketches of position- velocity diagrams of a blister flow HII region at three different 
viewing angles. The molecular cloud velocity is Yamh and the star is stationary. No stellar 
wind is present. The cloud material is continually evaporated by the stellar radiation and 
forced by pressures to flow along the density gradient. The Voo is the exiting speed of 
the ionized gas. The dashed lines show the locations of the p-v cuts. The figures show 
p-v diagrams of the region when the head of the region is tipped away from the observer 
{9 > 90°). The p-v diagrams would be mirrored about Vamb for stellar motions toward the 
observer (6 < 90°). 
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Table 1. OBSERVATIONAL DATA 



Object 


Date 


RA(J2000.0) 


DEC 


d 


^LSR 


SW 


SL 


STS 


PS 


ME 


IT 




mm/dd/yy 


hh mm ss.sss 


dd mm ss.sss 


(kpc) 


(km s-i) 


(") 


(") 


{") 


{") 


{" X ") 


(sec) 


2G5.89 -0.39 


06/11/03 


18 00 30.4 


-24 04 00 


2.6 =» 


9" 


0.9 


10.8 


0.4 


0.360 


6x6 


180 


iGll.94 -0.62 


07/02/03 


18 14 02.2 


-18 53 26 


4.2 " 


39 " 


1.4 


10.8 


0.7 


0.360 


18x16 


12 


2G29.96 -0.02 


06/11/01 


18 46 03.9 


-02 39 22 


7 A" 


98 <= 


0.9 


10.8 


0.4 


0.360 


20.4x13 


16 


^030. 54 +0.02 


06/30/03 


18 46 59.5 


-02 07 26 


13.8 '^ 


48 <= 


1.4 


10.8 


0.7 


0.360 


12x12 


16 


2G33.92 +0.11 


06/11/01 


18 52 50.2 


00 55 29 


8.3 " 


108 " 


0.9 


10.8 


0.4 


0.360 


17x15 


50 


2G43.89 -0.78 


07/01/03 


19 14 26.2 


09 22 34 


4.2 '' 


53.5 = 


1.4 


10.8 


0.7 


0.360 


12x12 


16 


2G45.07 +0.13 


07/01/03 


19 13 22.1 


10 50 53 


6.0 " 


59== 


1.4 


10.8 


0.7 


0.360 


7x8 


12 


2G45.12 +0.13 


07/01/03 


19 13 27.8 


10 53 36 


6.9 " 


59 " 


1.4 


10.8 


0.7 


0.360 


15x18 


16 


2G45.45 +0.06 


06/11/01 


19 14 21.4 


11 09 14 


6.6 == 


59.5 <= 


0.9 


10.8 


0.4 


0.360 


15x12 


52 


iW51 IRS2 


06/30/03 


19 23 39.9 


14 31 09 


6.6 <= 


58.3 " 


1.4 


10.8 


0.7 


0.360 


30x14 


12 


1G61.48 +0.09B 


07/01/03 


19 46 49.2 


25 12 43 


5.4 " 


22 = 


1.4 


10.8 


0.7 


0.360 


60x19 


16 


1K3-50A 


06/30/03 


20 01 45.7 


33 32 43 


8.7 f 


-24.4 « 


1.4 


10.8 


0.7 


0.360 


12x12 


24 


1S106 


06/29/03 


20 27 26.8 


37 22 48 


0.6 s 


-1.0 ^ 


1.4 


10.8 


0.7 


0.360 


47x53 


16 


1NGC7538A 


09/13/02 


23 13 45.6 


61 28 18 


3.5 ' 


-56.9 J 


1.4 


10.1 


1.0 


0.360 


29.4x16.2 


>20 


iW3A 


07/01/03 


02 25 40.8 


62 05 53 


2.0 "^ 


-40' 


1.4 


10.8 


0.7 


0.360 


86x37 


>8 


iW3B 


07/03/03 


02 25 36.9 


62 05 45 


2.0 ^ 


-43 ™ 


1.4 


10.8 


0.7 


0.360 


61x44 


>8 



Note. — The observational parameters (SW: slit width, SL: slit length, STS: step size, PS: plate scale, ME: map extent, IT: pixel integration 
time) of observations. The shown coordinates are for the (0,0) positions in [Ne II] line maps. Two methods are used in determining these 
coordinates: Read from published figures. Both qualities of the figures and the reading process can affect uncertainties of the resulting 
coordinates, which can be over 1". -^ Cross-correlate [Ne II] line maps resampled on fine grids (pixel size <0.1") and available radio continuum 
maps to find the best morphological fit. Uncertainties of these coordinates should be better than 0.5". Nasa/IPAC Extra gal actic Database 
NED ) c oordinate calculator is u sed to convert obtained B1950.0 co ordinates to J2000.0 coordinate s. '^_y<in^^<od | |2003[): ' 'iHatch ell et al.l 



1998h 
'' ISchneider et al 



Churchwell et al.l l|l990l) 

" ll2002h: ' iHanson et 



only for W3B, 



Wood fc Churchwelll lll989bl ' 

" l l2002h J iDickel et al.l Jl98ll 



Tieftrunk ct al.. 1199 



t al.l ( 

m 



IShepherd fc Churchwehl lll996l') : ^ iHarrisI lll9~ 

'^ iHachisuka et al.l l l2006h : ' only for W3A, iKantharia et al.l 



g IStaude et al.l 



1982) 



19981) : 
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Table 2. OBSERVED KINEMATIC TYPES 



Object 



Type Comments 



Mon R2 


C/D 


G5.89 -0.39 


7 


G11.94 -0.62 


7 


G29.96 -0.02 


D 


G30.54 +0.02 


C/D? 


G33.92 +0.11 


D 


G43.89 -0.78 


D 


G45.07 +0.13 


7 


G45.12 +0.13N 


7 


G45.12 +0.13SE 


7 


G45.12 +0.13SW 


7 


G45.45 +0.06 


7 


W51 IRS2 


D 


G61.48 +0.09B 


7 


K3-50A 


7 


S106 


7 


NGC7538A 


D 


W3A 


C/D 


W3B 


D? 



70° viewing angle, see ^hi^_et_aL (2005) 

barely resolved 

multiple peaks, extinction lanes, complex kinematics, V ~ Vamb 

140° viewing angle, sec Zhu ct al. (2005) 

horseshoe morphology, weak type D signature on axis 

confused morphology, but clear kinematic signature 

clear type D morphology and kinematics 

unresolved, 1" diameter 

appears cometary but without obvious cometary kinematics 

barely resolved 

broken ring, multiple peaks 

confused line morphology, with different [Ne II] and S morphologies 

[S IV] has e ^ 180°, type D kinematics (Lacy et al. 2007) . 

[No II] shows superimposed sources 

extinction lanes, complex line morphology and kinematics, V > Vmolec 

three blended peaks, extended line emission 

bipolar, with unclear kinematics 

shell-like appearance, tail-on type D kinematics 

cometary with 9 Ri 90° 

shcU-likc appearance, confusing kinematics 



Note. — Kinematic types: A. expanding shell, B. expanding sphere, C. bow shock from moving star, 
D. pressure-driven surface flow, E. blister. 



